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We present a set of notes, meant for quick reference, on radiative po- 
larization, computer algorithms and spin matching in electron storage 
rings. 



2.7.7. Radiative Polarization in Electron Storage Rings 

The Sokolov-Ternov effect [jl[] Relativistic electrons in a storage ring emit syn- 
chrotron radiation (Sec. 3.1 in [2~7|). A very small fraction of the radiated photons 



cause spin flip. For electron spins aligned along a uniform magnetic field, the f J. and 
It flip rates differ and this leads to a build-up of spin polarization antiparallel to 
the field. Positrons become polarized parallel to the field. The transition rates for 
electrons are 
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For positrons, interchange plus and minus signs here and elsewhere Q. 

The equilibrium polarization in a uniform magnetic field is independent of 7, 
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0.9238 (2) 



For a beam with zero initial polarization, the time dependence for build-up to equi- 
librium is 

Pit) = P st [1 - exp H/to)] (3) 

where the build-up rate is 
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a Contributions to the Handbook of Accelerator Physics and Engineering, Eds. A.W. Chao and M. 
Tigner, World Scientific, 1999. Also as DESY report 99-095. 

b p is the radius of curvature of the orbit and the other symbols have their usual meanings. 
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t depends strongly on 7 and p but is typically minutes or hours. When p depends 
on the distance along the design orbit s, i.e. p = p(s), just average Eq.(^) over the 
circumference C. Then for a flat ring in which all bending magnets have the same p 

[ J 99C[m]p[m]2 W 



The Baier-Katkov flip rate For electron spins initially aligned along an arbitrary 
unit vector £ the generalization of Eq. (U) is 



w = — 
2t„ 



(6) 



where s = direction of motion and b = (s x s)/|s|. b is the magnetic field direction if 
the electric field vanishes and the motion is perpendicular to the magnetic field. 
The corresponding instantaneous rate of build-up of polarization along £ is 



'bk 
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(7) 



The T-BMT equation Neglecting radiative spin flip, the motion of the centre-of- 
mass spin expectation value £ of a relativistic charged particle traveling in electric and 
magnetic fields is governed by the Thomas-BMT equation d^/dt = Q x £ (Sec. 2. 7.1 

in HP- 

We write 

Q = tt rc{ + uj imp + uj sh (8) 

where f2 ref is due to design fields on the design orbit and f2 ref +a? mp (= Q co ) is due to 
fields on the closed orbit whereby d? mp is due to field imperfections and corrections. 
uP b is due to synchrotron and/or betatron motion with respect to the closed orbit. 
On the closed orbit the T-BMT equation 

|f=rrxf ( 9 ) 

can be solved in the form 

^ = Rz 3 (s,s )a So ) (10) 

where R3X3 is a rotation matrix. The real unit eigenvector (rot. axis) for the one turn 
matrix R3° 3 (s + C, s), denoted by n Q (s), is the periodic spin solution on the closed 
orbit. For a perfectly aligned flat ring with no solenoids, n Q (s) = ±y. The one turn 
matrix has a complex conjugate pair of eigenvalues e ±l2nUspin . Given n , we introduce 
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a pair of unit vectors (m , Iq) such that n x m = / , / x n = m , fulfill Eq.(^), and 
obey 



m (s + C) + il {s + C) 



g227ri/ spin 



mo(so) + ^o(«o) 



The (mo, ^o) are usually not periodic in s. But by applying a further rotation by an 
angle ip sp \ n (s) around n we can construct the vectors (m, I), 



rh(s) + i l(s) 



-i^spin(s) 



m (s) +il (s) 



(12) 



By choosing ^ sp i n (s + C) — ip sp m(s) = 27rz/ spin , the set (h ,m, I) is then periodic in s 
with period C. The vectors (m, I) are needed in Sec 2.7.8. 

The closed orbit spin tune v sp \ n is the number of spin precessions per turn around 
n . For a perfectly aligned flat ring without solenoids z/ spin = 07, where a = (g — 2)/2 



(see Sec. 2. 7.1 in ]27]l). In this section and in Sec. 2.7.8 we use the symbol "a" instead 
of the symbol "G" used in the rest of the Handbook. Only the fractional part of the 
spin tune can be extracted from the numerical values of the eigenvalues e ±l2nl/spin . For 
the definition of spin tune away from the closed orbit see 0. 



The Baier-Katkov-Strakhovenko (BKS) equation Neglecting the effect of 
stochastic (synchrotron radiation) photon emission on the orbit and imagining that all 
particles remain on the closed orbit, the equation of motion for electron polarization 
is § 1 



dP 
~dt 



— = Q co x P 



(13) 



In the case of horizontal motion in a vertical magnetic field, we have Q = (a , yc/p)y, 
and b(s) = y. 

By integrating the BKS equation, one finds the generalized Sokolov-Ternov for- 
mula for the asymptotic electron polarization in arbitrary magnetic fields along the 
closed orbit, 



bks 
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See for a compilation of time scales. Usually, in rings containing dipole spin 
rotators (Sees. 2. 7.3, 2.7.4 in [pjj) the polarization |Pbks| cannot reach 0.9238 0. 
The BKS polarization build-up rate is 
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This is in accord with Eq.(^) by replacing £ — > ho and averaging. 
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Radiative depolarization The stochastic element of photon emission together 
with damping determines the equilibrium phase space density distribution. The same 
photon emission also imparts a stochastic element to u sh and then, via the T-BMT 
equation, spin diffusion (and thus depolarization) can occur 0. The polarization is 
the result of a balance between the Sokolov-Ternov effect and this radiative depo- 
larization. In the approximation that the orbital motion is linear, the value of the 
polarization is essentially the same at each point in phase space and azimuth and the 
polarization is aligned along the Derbenev-Kondratenko vector n M. 

The unit vector field h depends on s and u = (x,p x ,y,p y , z,5). h(u;s) satisfies 
the T-BMT equation at (u; s) and is periodic: n(u; s) = n(u; s + C). On the closed 
orbit n{u; s) reduces to n Q (s). 



The Derbenev-Kondratenko— Mane formula Taking into account radiative de- 
polarization due to photon-induced longitudinal recoils, the equilibrium electron po- 
larization along the h field is || [jl], [| 

P = 8 f ds (w i -<- fl -&,, 

where < > s denotes an average over phase space at azimuth s. This formula differs 
from Eq.(|l4|) by the inclusion of the terms with || and use of h instead of no- The 
ensemble average of the polarization is 

PeoB,ek(s) = Pdk (n) s (17) 

and (n) s is very nearly aligned along no(s) (see the angle estimate below). The value 
of the ensemble average, P e ns,dk(s), is essentially independent of s. 

The effect of transverse recoil can also be included but contributes derivative terms 
analogous to ^ which are typically a factor 7 smaller than ^ and can be neglected 
unless ^ is very small O, |T2|. If ^ were to vanish, a Pdk of 99.2 % could be reached 
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In the presence of radiative depolarization Eq.fllSD becomes 



r dk 



8 rale 2 C J \ \p(s)\ 
This can be written in the form: 



111 
— = - + , (19) 

7"dk Tst Tdep 



where r^ 1 can be (very well) approximated by r bks in (15) and 



r ^T^Mt> (20) 
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The time dependence for build-up from an initial polarization P to equilibrium is 

P{t) = Pen S ,dk [1 - e~^] + P e-* /Tdk (21) 

This formula can be used to calibrate polarimeters (see Eqs.(21) and (22), Sec. 2. 7. 8) 
fl"3|| . However, the calibration will be imprecise if ^ in the numerator of Eq.([Uj) is 
not well enough known. For examples of build-up curves see 0. 



Resonances Away from the spin-orbit resonances (see also Eq.(ll), Sec. 2. 7. 8) 

^spin = k -\- k x v x + kyfy + k z v z (22) 

n(u; s) ~ fio{s). But near resonances h(u; s) deviates from n (s) by typically tens of 
milliradians at a few tens of GeV and the deviation increases with distance in phase 
space from the closed orbit. The spin orbit coupling function ||, which quantifies the 
depolarization, can then be large and the equilibrium polarization can then be small. 
Note that even very close to resonances, \(n) s \ ~ 1: the ensemble average polarization 
is mainly influenced by the value of Pdk in Eq. ([16|) . 

To get high polarization, one must have (dn/d5) 2 <C 1 in dipole magnets. The 
machine optimization required to make || small is called spin matching (Sec. 2. 7.8). 



Asymmetric wigglers If t^ s is very low because the energy is low and/or the 
average curvature is small the polarization rate can be enhanced (see Eq.([L5|)) by 
installing an "asymmetric wiggler" , i.e. a string of dipoles in which short dipoles with 
high fields are interleaved with long dipoles with low fields of opposite polarity while 
ensuring that the field integral of the string vanishes. For more details, and discussion 
of advantages and disadvantages see 0. A particular potential disadvantage is that 
the enhanced radiation loss can require that extra rf power be installed and that the 
energy spread increases so that the depolarization rate increases owing to stronger 
synchrotron sideband resonances (Sec. 2. 7.8). 



Kinetic polarization The (numerator) term linear in ^ in Eq.(|T6|) is due to a 

correlation between the spin orientation and the radiation power ||. In rings where 
no is horizontal due, say, to the presence of a solenoid Siberian Snake (Sees. 2. 7. 3, 
2.7.4 in l27f) ||14|| , || has a vertical component in the dipole fields. This can lead to a 
build-up of polarization ( "kinetic polarization") even though the pure Sokolov-Ternov 



effect vanishes. 



The rate is r^ 1 . 



Phase space and polarization evolution equations If the orbital phase space 
density ip obeys an equation of the Fokker-Planck type (Sec. 2. 5. 4 in p7| ) 

f = (23) 
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where £ pp is the orbital Fokker-Planck operator, then the spin diffusion is described 
by the "Bloch" equation 

^ = £ FP V + nxV (24) 
as 

where fl = Q/(ds/dt) and V is the "polarization density" = 2/hx (density in phase 



space per particle of spin angular momentum) [15|, 16]. To include the Sokolov-Ternov 
effect see IffTI. 



Beam energy calibration A polarized electron beam can be depolarized by apply- 
ing a weak oscillating magnetic field perpendicular to n with a frequency f r f related 
to the fractional part of the spin tune z> sp i n by 

/rf = /c^spin Or f ri = / c (l - Z>spin) (25) 



where f c is the circulation frequency of the beam [|18[]. Thus the required f T { gives an 
accurate measurement of z> sp i n and this gives high relative precision knowledge of f sp i n . 
By relating z/ spin to the average energy of each beam, high precision measurements 
of the centre-of-mass energy of colliding e + -e~ beams and of the masses of vector 



mesons such as the T's and the Z can be obtained [T^, ^T], ^2| |23| . Other beam 
parameters can also be measured pi"| . The polarization need not be large for these 
measurements so that by Eq.(21) the depolarization can be repeated at intervals of 
about Tdk- 



Concluding remarks For an overview of measurements see [p5j, |7], (26|. For an 
overview of the theoretical background see |3| . 
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2.7.8. Computer Algorithms and Spin Matching 

There are two classes of computer algorithm for estimating the equilibrium polariza- 
tion in real rings: 

(i) Methods based on evaluating || in the Derbenev-Kondratenko-Mane (DKM) 
formula (Eq.(16) of Sec. 2. 7. 7) given the ring layout and magnet strengths; and 

(ii) The SITROS [I] algorithm which estimates Td ep (Sec. 2. 7. 7) using Monte-Carlo 
tracking. 

The class (i) algorithms are further divided according to the degree of linearization 
of the spin and orbital motion: 

(ia) The SLIM family (SLIM |, |, SLICK [|, SITF 0) and SOM § and ASPIR- 
RIN ||. The latter two utilize the "betatron-dispersion" formalism outlined below 
and all are based on a linearization of the orbital and spin motion. 

(ib) SMILE [0]: Linearized orbital motion but nonlinear spin motion; 

(ic) SODOM ||: Linearized orbital motion but nonlinear spin motion; 

(id) SpinLie: Nonlinear orbital motion and nonlinear spin motion (Sec. 2. 7. 9 in 
H); and 

(if) SPRINT ||: Linearized orbital motion but nonlinear spin motion. 



The linear approximation — SLIM We now present expressions for ^ in an 
approximation in which the orbit and spin motion are linearized and in which u; sb 
(Sec. 2. 7. 7) is linearized as in Eq.(2) below (the SLIM formalism). In linear approxi- 
mation we write (see Sec. 2. 7. 7) 

h(u; s) = n (s) + a(u; s)rh(s) + (3{u; s)l(s) (1) 



valid for \J 'a 2 + (3 2 <C 1 and we write the components oof*, u>^, in the form 



ltt 11 



cof> 



3x6 



UJ. 
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/ X \ 

Px 
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Py 

z 

\ & ) 



(2) 



where u = (x,p x ,y,p y , z,5) describes motion with respect to the closed orbit. In 
particular p x = x' and p y = y' (except in solenoids) . 

The detailed forms of the matrix F 3x6 for bending magnets, quadrupoles, skew 

The orbit motion 



quadrupoles, solenoids and rf cavities can be found in [11, 12 



in sextupoles is linearized. For example for a quadrupole, defining g 
where q = one has 

•> po ox 



•(i + ai)g 



F(s) 




0^000 
(? 



(3) 
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In linear approximation the combined orbit and spin motion is described by 8 x 
8 transport matrices of the form 

M = ( ™ 6x6 ° 6x2 ) (4) 

acting on the vector (u, a, (3), where M 6X 6 is a symplectic matrix describing orbital 
motion and G2X6 describes the coupling of the spin variables (a, (3) to the orbit and 
depends on m(s) and l(s) (see e.g. Eq.(14)). D 2X 2 is a rotation matrix associated 



with the spin basis rotation of Eq.(12) in Sec. 2. 7. 7 [[H], PH . 

The eigenvectors for one turn defined by M(s + C, s ) ■ q^ = A M • g M are written 
in the form 

Qk(s ) = (j ( ( S ; o ))> *-*(*o) = [&(*>)]• 
for k = I, II, III ; 



6 (s ) 



, q-k(so) 



[qk(so)Y 



for k — IV 



(5) 



for arbitrary s - The Vk are the eigenvectors for orbital motion with eigenvalues 
Afc = e~~ i^nvk anc j v _ k _ _ Vk ^ _ j ? jj^ j/j). These eigenvectors obey 

the orthogonality relations, and have the normalization of 0. The corresponding 
eigenvalues of M(s + C, s ) are = A^ (k = I, II, III) and A JV = e — with 



v IV = z^spin and with u_ IV 



-is, 



The spin parts of the eigenvectors Wk(so) (k = I, II, III) and w iv (sq) can be 
written as 



Wk{S ) 



D(s + C, s ) - X k G(s + C,s )v k (s ) 



for k — I, II, III ; 



(6) 



w IV {s ) 



_L_ / 1 

V2 



-»^spin(S0) 



for fc = IV 



and 



w„ k (s ) = [w k (s ))*,(k = I,II,III,IV) 



In this linear approximation h{u; s) can be obtained via |13|, [14|, |TT 

h(u; s) - h (s) = ^ ^ J = + A_ fc t2L fc (s)} 

9 



(7) 



where the amplitudes A k are determined by the orbit via 

(s) = {A k Vk(s) + A-kV-kis)} (8) 



u 



Then with respect to the (h ,rh,l) frame, 

dh 



85 



= i Yl Ms^k ~ v k5 w* k } 

k=I,II,III 

= -2 1m v ^ (9) 

k=I,II,III 

Note that this is independent of the phase space vector u and that || is periodic in 
azimuth in the machine coordinate system. In this approximation the depolarization 
time is then (Eq.(20), Sec.2.7.7) 



_! 55>/3e 2 7 5 /il r o+c ,_ 1 r * ,~\ 



(10) 



This is the formula used in SLIM to calculate the depolarization rate. SLIM is 
based on thin lens optics. SLIM-like programs for thick lens optics are SLICK and 
SITF. 

In this linearized theory the vectors n(u; s) and || display only first order reso- 
nance behaviour, namely the resonances 

^spin = k + kiui + k n v u + kiiiViu (11) 

with \kj\ + \kn\ + \ kjjj\ = 1. They arise from the denominator matrix in Eq.(6). The 
theory is not valid beyond the limit a/q; 2 + /3 2 < 1. In a fully uncoupled optic I — > x, 
II — > y and III — > z. 

This formalism forms the natural language for the method of maximizing the 
polarization called "spin matching". Thus comments on the other programs will be 
postponed until later. 



Spin matching in the SLIM formalism In practice the spin matching of real 
rings takes place in stages as follows. 

Stage 1: Strong synchrobeta spin matching of the perfectly aligned ring 

From Sec.2.7.7 it is clear that to maximize the polarization we must minimize r^. 
Then by Eq.(10) we need to minimize v kb (k = I, II, III) or the components of w k 
at azimuths where l/|p(s)| 3 is large. The vl 5 determine the orbit excitation due to 
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synchrotron radiation (Sec. 2. 1.4 in [[HJ. In particular, for rings without x — y 

coupling, v* jr usually vanishes in the arcs since the vertical dispersion r\ n vanishes. 
However, v* jr does not vanish inside spin rotators (Sees. 2. 7.3, 2.7.4 in ||45|| ) contain- 
ing vertical bends. On the other hand v* tends not to vanish in the arcs since the 
horizontal dispersion n 1 ^ 0. Finally, u* essentially never vanishes. Each case 
must be evaluated individually but the minimal recipe is to try to minimize Wk for 
{k = I, II, III) only at azimuths where |f fc5 (s)| 2 /|p(s)| 3 is sufficiently large. This 
in turn requires (Eq.(6)) that G(s + C, s) ■ Vk( s ) f° r (k — !■> H, HI) be minimized. 
This must be achieved by designing the ring layout with this in mind and then pro- 
viding sufficient flexibility in the optics by providing enough independently powered 
quadrupoles. Subsequent calculations with SLIM will indicate whether the match 
criteria for the adopted design suffice. 

Consider, for example, a specific mode, k. Label those bending magnets at which 
|f fc5 (s)| 2 /|p(s)| 3 is large by Hi^ k \ //2 > • • l^n^ ■ Then the suppression of depolar- 
ization associated with the kth mode requires that w/^s^J = for all (i — 1 to n*.). 



In general (see Eq.(6)) this in turn requires |T6| 



G(s M2 , s^)v k {s^) = 
G(s M , s w ) v fc (s w ) = 

G^+Cs^v^) = (12) 

where we suppressed the superscript label u k" . To fulfill Eq.(12) we then require the 
Gij(sn l+1 , s M; ) to vanish when the jth component of Vk does not vanish. The matrix 
G can be written in the form 

G(s 2 , si) = / dsT>(s 2 , s)G (s)M(s, si) 

J Si 



where 



Go = ( h lx ly ) • F (13) 



-m s -m x -m y 



Thus Gij(s m+11 s M; ) depends on the orientation of the (m, /) vectors so that in some 
cases some elements of Gij(s N+1 , s w ) vanish automatically. But in general these con- 
ditions can only be fulfilled by adjusting quadrupole strengths — while maintaining 
other necessary features of the orbital optics. We call this strong synchrobeta spin 
matching. A section of the ring satisfying a condition in Eq.(12) is "spin transpar- 
ent" for mode k. The interpretation is immediate: the overall spin-orbit coupling for 
the section vanishes for mode k. Clearly, the exact spin matching conditions are very 
dependent on the layout of a machine and each case must be handled individually. 
In thin lens approximation the G matrix for a quadrupole of length l q is 

/ -q ty o -ql x \ 

\ +qm y +qm x / 1 ' 
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where q = (1 + cry) g l q - The thin and thick lens forms of G for other magnet types 
are given in |3|, |17|] . 

If the Gij(s N+1 , s w ) cannot be brought to zero while maintaining an acceptable 
optic, then the G(s w + C, s w ) ■ Ufc(s w ) themselves should be minimized. This essen- 
tially means that the effects of elements of the G matrices of sections of the ring are 
made to partially cancel one another. The spin matching of a ring with a solenoid 
Siberian Snake (Sees. 2. 7.3, 2.7.4 in f4"5|l ) has provided an example of this 0. By 
Eq.(7) reduction of G(s + C, s) ■ v k (s) f° r {k = I, II , HI) also reduces the angle 
between n and n at azimuth s. 

Alternative Stage 1: Harmonic synchrobeta spin matching of the perfectly aligned ring 

If the strong spin matching methods just described are impractical for some rea- 
son, another approach aimed at minimizing the strengths of depolarizing resonances 
can be adopted. 

Rewrite Eq.(6) as 



3 ±#s P in(so+C) rso+C 



ds-j^\§)e- i27T l Uk±Us ^Y s/c 



so 



with 

j^\s) = e ±i[2 ™^~ s/c -^ m (l s ±im s l x ±im x l y ±im y ) Fv k (s)e +i27TUk ' s/c 

+00 

p=— 00 



c (T) - 
fcp ^ 



X - [ d§e i2w i Uk±Us ^- p Y /c e Tii, ^ n{ ' s) ( Z s ±im s l x ±im x l y ±im y ) F4(5) 
^ Jo 



so that 



f» +°° -i27r[i' fc ±i/ S pi n -p]so/C 

w kl (s )Tiw k2 (s )} = e ^^o) J- eg 



p=-oo ^ - ' s P ln 



The condition that w k (s N ) = for all (z = 1 to n k ) is now be replaced by 

r* +°° „-«27r[^±i/ S p in -pls^./C 



2vr kv [u k ±u spin -p] 

Near to the resonance v k ± z/ spin — p = the sum over p is dominated by the 
term containing This corresponds to the spins' seeing a stationary field in the 
(ho, 777.0) ^0) frame, proportional to which rotates spins away from ho- Note that 
c^p is independent of s^. 
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Approximate spin matching can be achieved for all s Mi by adjusting the optics so 
that an appropriate set of the are small. This is called harmonic synchrobeta 
spin matching. See also |L6], Q. For mode k and orbit amplitudes A k and A- k 
the strengths of the resonances are proportional to Af.c^~ and A_fcC^~. The concept 
of resonance strength is important for the acceleration of polarized protons (Sec. 2. 7. 5 
in Normally only the case of flat rings is considered so that no is vertical. 

The formalism presented here shows how to define resonance strength for arbitrary 
orientations of n . See also [0| |20|. On resonance with e l2 A Vk±v ^ m ~p] s / c = 1, the 
coefficients take the form 



1 

C 



/ dse" #spin(s) ( l s + im s l x + im x l y + im y ) • Fv k (s) 
Jo 



for v k + z/< 



spin 



P 



C (+) 



1 r c 

- / dse +i ^ { ~ s) ( l s - im s l x - im x l y - im y ) • Fv k (s) 
^ Jo 



C 

for v k - i/ spin 



V 



(15) 



These can be obtained from the SLIM algorithm by calculating the matrix G at 
the resonance for one turn but without the backward spin basis rotation (Eq.(12), 
Sec. 2. 7. 7) that, in SLIM, is applied at the end of one turn 



Reformulation in terms of beta functions and dispersion^A\ 



We can reformulate Stage 1 by making a transformation of the particle coordinates 
from u = (x, p x , y, p y , z, 6) to u = (x , p x , y, p y , z, 5) via the transformation 



u 



K-u 



where 



K(s) 



( 1 





\ 





1 
1 


-r)i ^4 









1 










~V3 1 





-V3 



1 / 



whereby the dispersion vector fj = (rji, r] 2 , r] 3 , 774) is the periodic solution of the lin- 
earized equations of motion for (x,p x ,y,p y ) with 5 = 1 and without the rf cavities. 
Then 



x 



x-Srjx, y = y-Sr/2. 
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The matrix K is symplectic so that the formalism remains canonical. In particular, 
the new transfer matrices M and eigenvectors are obtained via 

M(s 2 ,si) = K(s 2 ).M(s 2) si).K-V) 

and 

M(s + C,s) = K(s) ■ M(s + C, s) ■ K _1 (s) 



K( S )^( S ) 



so that the eigenvalues and orthogonality conditions are unchanged fl2| . Furthermore 
the new matrices F and G are 



F(s) = F(s)-K" 1 (s) 

and 

G(s 2 , Sl ) = GKsO-K- 1 ^) 
The depolarization rate then takes the form 

with /fc = 5^n=i [-^ _1 ]5n ' ^fcn = ^5 an d ^fe = Wk- This formulation has the advantage 
that in the special case, or the approximation, of no orbital coupling, the 6x6 orbit 
matrices just consist of three 2x2 matrices on the diagonal. 

The eigenvectors Vk{s) of the revolution matrix are then written in the form 




(r = x,y,z) and the f k are given by fj = f x = -{vrin* ~ Vrinx) 5 /// = fy = 
-(v II3 T]4 - VhaVs) and fm(s) = f z = J~^-e~ l ^ s \ The functions // and f n are 
closely related to the factor 

Vk 2 + {am + PkV' k ) 2 ( k = j jj) 
14 



in P2j. f in is in practice almost independent of s (see below). Note that these 



a's and /3's are Courant-Snyder parameters and should not be confused with the 
quantities in Eq.(l). With these coordinates the F matrix for a quadrupole takes the 
form 


g £773 
g g Vl 

We can write 

-l _ 



w k (s 



GOo + C, s ) ■ v k (s ) 



- - D(s + C, s ) - X k 
for (k = I, II, III) and we use a representation of the G matrix in the form 
G(s 2 , si) = / D(s 2 , s) G (s) M(s, si) 



with 



G„ [ ls lx ly 1 ■ F 

-m, — m T — m„ 



In thin lens approximation the G matrix for a quadrupole is 

/ -gij, «i 

V +qm y +gm x k 2 

where K\ = —ql y rji — ql x r] 3 and k 2 = +qm y rji + qm x rj 3 . We see that as a result 
of separating the transverse coordinates into betatron and dispersion contributions, 
columns six of F and G contain terms depending on dispersions. 

The strong spin matching condition Wk — for suppressing depolarization now 
amounts to setting the G(s Mi , s rt ) uj.(s tt ) to zero in analogy with Eq.(12). Then 
in the special case, or approximation, of a fully uncoupled optic and by taking into 
account only the depolarizing influence of quadrupoles this is equivalent to requiring 

mm-- 



For horizontal motion: 

_ (1 + ot) 1 p +1 rf5 y^^ (5)e -^4 S -)[ /v(s ~ )± , mv(s ~ )]e T^ Pin « = o 

V 2 c 



For vertical motion: 

(1 + 07) 1 /" s ^+i 

71 - c 



(16) 



dS^P y ~{5)g(s)e~ 1 ^ [l x (s) ±im x (s) } e ^-«« = 

(17) 
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For longitudinal motion: 



;i + Q7) i 

y/2 C 



,Ja z {S)+i\ 
as 



-9{s)e 



x{ Vy [l x (s) ± im x (s)\ + Vx [l y (s) ± im y (s)]}e Tl ^ {s) = (18) 

Since in practice synchrotron motion is well approximated by simple harmonic 
motion [25|, f3 z (s) is almost independent of s and ct z (s) ~ 0. Then Eq.(18) may be 
approximated by 



1 + 07) i i r^+i 



dsg(s)e 



-iipz(s) 



X 



V2 VW z c 

; T^ P in(s) {% ±i mx (s)} + r] x [l y {s) ±im y (s)}} 







(19) 



Harmonic synchrobeta spin matching in terms of beta functions and dispersion 
follows the path detailed earlier under "Alternative Stage 1" but with the eigenvectors 



Vk and the matrices F. Typical expressions can be found in [24 



Commentary 



Spin matching should be carried out using thick lenses so that the optic is correct. 
Strong spin matching by minimizing the integrals in Eqs. (16-18) requires explicit in- 
tegration. Furthermore Eqs. (16-18) must be modified if there is significant orbital 
coupling. Thus in practice the numerical fitting involved in strong spin matching 
can be carried out most simply by minimizing the GV,(s Mi+1 , s M; ) since these already 
represent integrals and do not need knowledge of the dispersion. Moreover these 
matrices are precisely those contained in the SLIM program so that cross checks be- 
tween programs are simplified. Another advantage of working with the G matrix is 
that it allows sections of the ring to be studied and made transparent in isolation 
since no knowledge of Courant-Snyder parameters is needed; use of G emphasizes 
the local nature of spin transparency. On the other hand Eqs. (16-18) and the split- 
up versions depend on Courant-Snyder parameters and these in turn depend on the 
structure of the whole ring so that the "locality" is masked. When studying the spin 
transparency of a ring, it is often useful for diagnostic purposes to set elements of the 
G or the G matrices to zero artificially and thereby obtain an impression of which 
sections of the ring are most dangerous. For example by switching off column six of 
G in quadrupoles, the effect of dispersion can be cleanly separated from the effect of 
betatron motion. One can also investigate the system by using the matrix handling 
facilities in symbolic algebra programs and the fact that the G and G of magnets or 
strings of magnets often depend in a simple way on the elements of the corresponding 
M and M [fT/]] . Finally, the G and G matrices are in general energy dependent. But 
a spin match made at the design energy is usually still effective for a few tens of MeV 
above and below, except near resonances. 
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Some examples 



In a perfectly aligned flat ring (no vertical bends) with no solenoids and no x — y 
coupling, the depolarization rate Un vanishes (see below) so that no spin matching 
is needed. 

A spin rotator (Sees. 2. 7.3, 2.7.4 in JIIJ) based on dipoles and containing no 
quadrupoles is automatically almost spin transparent since the elements of G are usu- 
ally smaller in dipoles than in quadrupoles [[H]]. Dipole rotators containing quadrupoles 



need explicit spin matching 

Spin rotators based on a combination of solenoids (which rotate no from the 
vertical into the horizontal) and dipoles (to make the polarization longitudinal at an 
IP) (T7J are not automatically transparent. They also cause x — y coupling. However, 
by sandwiching quadrupoles and skew quadrupoles among sections of solenoid the 
coupling can be eliminated and by careful choice of the sandwich structure some 
terms in columns 1 to 4 of G for the rotator can be made small at the same time 



17fl . Column 6 remains troublesome but for antisymmetric solenoid schemes |L7| the 



columns 6 of the rotators cancel each other. For further discussion on solenoids see 



27, 28 



For a straight section (e.g. surrounding an IP) where the polarization is longi- 
tudinal and which only contains quadrupoles and drifts, the spin precession angle is 
a linear combination of the overall orbit deflections Ap x and Ap y in the quadrupole 
fields [ j~7| . Thus spin transparency implies making Ap x and Ap y vanish for all or- 



bits. If the straight section is geometrically and optically left-right symmetric, this 
can be achieved with an optic for which tan.Aip x = — a x and td,nAip y = — a y where 
the A-?/>'s are the phase advances between the IP and the outer end of the straight 
section and the a's appertain to the outer end. So the eight conditions that columns 
1 to 4 of the G matrix vanish have been reduced to two conditions by the symmetry. 
Furthermore, this is an example where the spin matching conditions reduce to purely 
optical conditions. 

These conditions can also be formulated directly in terms of G. By choosing 
I = y and m = x and requiring that the elements G\\ and G23 vanish for the stretch 
from the IP to the outer end, G vanishes for the whole straight section for arbitrary 
orientation of m, I around the longitudinal no- 

If the straight section contains rf cavities, their influence on the spin transparency 
can often be neglected. 

Other examples of the use of symmetry to simplify the spin match can be found 
in |23] where spin matching using variants of Eqs. (16-18) for a ring with dipole ro- 
tators is discussed. The results of a calculation with SLICK before and after a spin 
match can be found in |[29|| . Experimental observations resulting from successful spin 



matching involving spin rotators are described in pO]. 
Computer programs for strong spin matching 

Strong spin matching facilities based on evaluation of spin-orbit integrals (e.g. Eqs. (16- 
18) ) are built into the programs ASPIRRIN and SOM. To do spin matching in terms 
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of G the code SPINOR Bill can be used. 



Stage 2: Harmonic closed orbit spin matching 

Once the perfectly aligned ring has been spin matched, the effects of misalignment 
must be addressed. In a perfectly aligned flat ring with no solenoids, n is vertical so 
that l y and m y are zero. Then by inspection of the G matrix elements for horizontal 
bends, quadrupoles and rf cavities it is clear that for no x — y coupling, columns 1, 
2, 5 and 6 of G(s + C, s) vanish so that by Eq.(6) wi(s) and wiu(s) are zero. In 
particular columns 1 and 2 of G and columns 1, 2 and 6 of G vanish. Moreover, for 
no x — y coupling v} 15 vanishes. Then by Eq.(lO) T^ e lin is automatically zero. In 
rings with vertical bends (e.g. in spin rotators) uq is made vertical in the arcs by 
design. 

In real misaligned rings there is a vertical closed orbit distortion and n is tilted 
from the vertical in the arcs (see below) so that the above mentioned columns of 
G and G for the arc quadrupoles do not vanish. In practice the tilts can be tens 
of milliradians and they increase with energy (they are roughly proportional to try) 
but even these small angles can lead to strong depolarization so that it is essential 
that the ring be very well aligned from the beginning. Note that vertical closed orbit 
distortion leads primarily to depolarization due to horizontal synchrobetatron motion 
in the arcs. Note also that tilts of tens of milliradians cause a negligible decrease of 
the underlying ST polarization (Eq.(14), Sec. 2. 7. 7). 

If there is a vertical correction coil and a BPM near each quadrupole, one can 
try to minimize the combined vertical kick ( "kick minimization" ) [[32|] applied to the 
orbit by each quadrupole and its correction coil and thereby reduce the tilt of no due 
to the distorted orbit's being off centre in the (misaligned) quadrupoles. This also 
reduces the generation of spurious vertical dispersion so that the driving of v y and 
v z (Eq.(ll)) resonances is avoided. This presupposes that the positions with respect 
to the quadrupoles of the BPM's are well known. These relative positions can be 



estimated using beam-based calibration (Sec. 4. 5. 5 in [^5|)[p2|. However, kick mini- 
mization will not be effective if, say, the dipoles have significant tilt misalignments. 

If these measures are insufficient, a further method for bringing ho closer to the 
vertical is needed. n , and thus its tilt, for the distorted ring can be obtained as 
described in Sec. 2. 7. 7 but one gains more insight by using a perturbation theory 
based on SLIM concepts [Q. Viewed from the (h , m, I) frame calculated for the 
design orbit, the first order deviation of h Q from the design orientation can be written 
as 

[8n m {s) - i5n Q2 (s)] = -i— V, h k - 



'2ix^ k-u. 



k 



spin 



where the hf. are Fourier coefficients given by 

so+C 



1 rsa-t<^ 

h k = - I dS[di(S) - id 2 (s)} e- lk2 ™/° 



so 



Here 



a B 1+370 

di \ / I I I \ I / 




A J B,(l + a 7o ) 

where the AB XjVjS are field errors and w co is the deviation of the 6-D closed orbit from 
the design orbit. 5fiQ can be minimized by using correction coils to adjust the closed 
orbit (e.g. by generating closed bumps so that the luminosity is not affected) in such 
a way that the real and imaginary parts of hk, with k near z/ sp m, are small. This 
technique is called harmonic closed orbit spin matching and is embodied in the 



program FIDO |34|, |35[. See |L8[ also. If the machine distortions are not well known 
and if the closed orbit cannot be measured well enough, the closed orbit correction 
must be carried out empirically by observing the polarization. If the distortions 
and the orbit are well enough known the correction coil strength can be calculated 
ab initio (deterministic harmonic closed orbit spin matching) ||36|| . The correction 
scheme should be chosen so that it achieves the maximum effect on 5ho with the 
smallest possible additional orbit distortion. 

Harmonic closed orbit spin matching can in principle be used to minimize the 
5fiQ due to an uncompensated solenoid placed at the position of a nominally vertical 
no- However, this is achieved more efficiently by generating relatively antisymmetric 
vertical orbit bumps (spanning horizontal bend magnets) on each side of the solenoid 
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It might also be useful to weight 5n (s) by a periodic function p(s) In that 
case one tries to minimize p(s)5h (s). This is worth trying, for example, if the main 
source of depolarization due to misalignments is the coupling of non-zero l y and m y 
to the horizontal dispersion in the arcs (see Eq.(19)). This is often the case, as can be 
seen by examining the numerical values of the contributions of each mode (/, //, 
in Eq.(10). Then p(s) is taken to be rji. 

To minimize p(s)5n (s) one must minimize the harmonics hk of 

h(s) = p(s)(d l - id 2 ) + p'(s)[5n m (s) -idn 02 (s)] 
= ~h{s + C) 



whereby 



p(s) [6n 01 (s) - iSn 02 (s)} = Efc h 



k- 



Stage 3: Further tuning 

Harmonic closed orbit spin matching can generate spurious vertical dispersion and 
this in turn generates vertical emittance (nonzero v* IIb (Sec. 2. 1.4 in |[45|)) and also 
ensures that column 6 of G for the quadrupoles does not vanish. Thus extra depo- 
larization can occur. It might then be useful to overlay a harmonic vertical betatron 
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match (k = II in Eq.(15)) on any existing Stage 1 match, assuming that is possi- 
ble. Likewise, to overcome the effect of spurious vertical dispersion in column 6 of 
G one could use extra vertical correction coils to overlay a harmonic vertical disper- 
sion match (k = III in Eq.(15)). Usually both of these two extra matches would be 
empirical. One could also try to combine the harmonic closed orbit match and the 
harmonic vertical dispersion match into one procedure. 

Stage 4: Beam-beam spin matching 

The beam-beam interaction is equivalent to a nonlinear lens and can spoil a spin 
match. The effect of the beam-beam interaction on the polarization is not fully 
understood but it has been suggested that the beam-beam depolarization can be 
reduced by balancing the beam-beam deflection of spins against subsequent deflec- 
tions taking place in the ring quadrupoles. The condition for minimizing the effect of 
vertical kicks is independent of the current and charge distribution in the opposing 
beam and takes the form 



V Hy _j_ ^ Sill 2 ^ 



An equivalent prescription in SLIM formalism allows an arbitrarily coupled optic to 
be treated pl| . 



Higher order resonances To go beyond the linearization of spin contained in 
Eq.(l) one writes 

n{u;s) = (1 -a 2 - (3 2 ) 1/2 h (s) + am(s) + pt(s) (20) 

(for a 2 + f3 2 < 1) and does not linearize the T-BMT equation. Then spin-orbit 
resonances of arbitrarily high order can appear in ^ @. The strength decreases 
with the order (= \kj\ + \ku\ + \kjn\). In practice the most intrusive higher order 
resonances are those for which t/ spill = k ± i>\. + kjjjUjjj. These "synchrotron sideband 
resonances" of the first order parent resonances are due to modulation by energy 
oscillations of the instantaneous rate of spin precession around no- They originate in 
the part due to synchrotron motion in the term o> sb -n appearing in the full equations 
of spin motion (i.e. beyond the SLIM level) |4!|]. The depolarization rate associated 
with sidebands of isolated parent resonances (z/ spin = k ± v^) is approximately pro- 
portional to the depolarization rate for the parent resonances. Thus the effects of 
synchrotron sideband resonances can be reduced by doing the spin matches described 
above. Explicit formulae for the proportionality constants {"enhancement factors") 



can be found in [43, 44]. The underlying strength parameter (the "modulation index") 
of synchrotron sideband resonances is (ajas/uz) 2 which increases strongly with the 
energy and energy spread. 
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Other computer codes The SMILE algorithm is restricted to linearized orbital 
motion and calculates || by an extension of the first order perturbation theory of 
SLIM to high order using Eq.(20) and full 3-D spin motion. The algorithm involves 
multi-turn spin-orbit tracking. High order resonance effects are manifested by reso- 
nance denominators but the formalism ensures that the vector h is of unit length. The 
highest required absolute values of the kj, kjj, k HI are specified as input parameters. 

SODOM represents n by a spinor notation. The periodicity condition n(u; s) = 
n(u;s + C) (Sec. 2. 7. 7) is equivalent to periodicity in the three phases of linearized 
orbital motion and the one turn 2x2 spinor transfer matrix on a synchrobeta orbit is 
also periodic in the initial orbital phases. The spinor transfer matrix and n(u\ s) are 
then represented by Fourier series. The Fourier coefficients are obtained numerically 
and n(u\ s) can then be reconstructed. By constructing n at many points in phase 
space || can be obtained by numerical differentiation. The highest required absolute 
values of the ki, ku, km are specified as input parameters. 

The algorithm SpinLie utilizes Lie algebraic methods (Sec. 2. 7. 9 in |4"5| ) and can 
handle non-linear orbit motion and 3-D spin motion. 

The vector n(u; s) can also be obtained by stroboscopic averaging using the code 
SPRINT. ^ can then be calculated by numerical differentiation. This algorithm 
automatically includes all orders of resonance. 

The above algorithms all exploit the DKM formula (Eq.(16),Sec.2.7.7) but the 
SITROS program simulates the depolarization process directly using a Monte-Carlo 
tracking simulation of the effects on the orbit and then on the spin of stochastic photon 
emission and damping and delivers an estimate of Td op . The equilibrium polarization 
is then obtained from the approximation (Sec. 2. 7. 7) 

-Peq = -Pbks-^ 1 (21) 
Tbks 

where 

± = J- + -L. (22) 

r tot ^bks r dep 

This ignores the (normally small) term || in the numerator of the DKM formula. 



References 

1. J. Kewisch et al, Phys.Rev.Lett. 62 N°4 (1989) 419; The program SITF is part 
of the SITROS package. 

2. A.W. Chao, Nucl.Inst.Meth. 180 (1981) 29 

3. A.W. Chao, AIP Proc. 87, (1981) 395 

4. SLICK is a thick lens version of SLIM by D.P. Barber. Private notes (1982) 

5. K. Yokoya, User's manual of SOM: Spin-Orbit Matching (1996) 

6. C.W. de Jager, V. Ptitsin, Yu.M. Shatunov, Proc. 12th Int. Symp. High 
Energy Spin Physics, World Scientific (1997) 

7. S.R. Mane, Phys.Rev. A36 (1987) 105-130 



21 



8. K. Yokoya, KEK Report 92-6 (1992) 

K. Yokoya, DESY Report 99-006 (1999) and at the Los Alamos archive: 
physics/9902068 

9. K. Heinemann, G.H. Hoffstatter, Phys. Rev. E54 N°4 (1996) 4240 

10. H. Mais, G. Ripken, DESY Report 83-62 (1983) 

11. D.P. Barber, K. Heinemann, G. Ripken, Z. f. Physik C64 (1994) 117-167 

12. D.P. Barber, G. Ripken, DESY report on spin dynamics and spin matching in 
preparation. 

13. S.R. Mane, DESY Report 85-125 (1985) and Cornell Report CLNS-86/750 

14. D.P. Barber, K. Heinemann, G. Ripken, DESY Report M-92-04 (1992) 

15. A.W. Chao, J. Appl. Phys. 50 (1975) 592 

16. H. Mais, G. Ripken, DESY Report M-84-04 (1984) 

17. D.P. Barber et al, Particle Accelerators 17 (1985) 243 

18. R. Rossmanith, R. Schmidt, Nucl.Inst.Meth. A236 (1985) 231 

19. G.H. Hoffstatter et al., in Nonlinear and Stochastic Beam Dynamics in Accel- 
erators, DESY Report 97-161 (1997) 

20. V. Balandin, N. Golubeva, D.P. Barber, DESY Report M-96-04 (1996) 

21. H. Mais, G. Ripken, DESY Report 86-29 (1986) 

22. M. Sands, SLAC 121 (1970) 

23. A.W. Chao, K. Yokoya, KEK Report 81-7 (1981) 

24. J. Buon, K. Steffen, Nucl.Inst.Meth. A245 (1986) 248 

25. D.P. Barber et al, DESY Report 86-147 (1986) 

26. A. Blondel, J. Jowett, Proc. 8th Int. Symp. High Energy Spin Physics, AIP 
Proc. 187 (1989) 

27. A. A. Zholents, V.N. Litvinenko, BINP (Novosibirsk) Preprint 81-80 (1981). 
English translation: DESY Report L-Trans 289 (1984) 

28. D.P. Barber et al, DESY Report 82-76 (1982) 

29. D.P. Barber, Proc. 9th Int. Symp. High Energy Spin Physics, Springer (1991) 

30. D.P. Barber et al, Phys.Lett. 343B (1995) 436 

31. L. Hand and A. Skuja, DESY HERA Report 85-18 (1985); The program 
SPINOR has been updated and modernized by D.P. Barber. 

32. D.P. Barber et al., European Part.Acc.Conf. 1996 

33. D.P. Barber et al, DESY Report 85-44 (1985) 

34. S.R. Mane, DESY HERA Report 85-23 (1985) 

35. D.P. Barber et al, Nucl.Inst.Meth. A338 (1994) 166 

36. R. Assmann et al., Proc. 11th Int. Symp. High Energy Spin Physics, AIP 
Proc. 343 (1995) 

37. K. Steffen, DESY Report M-82-25 (1982) 

38. A. Blondel, CERN LEP Note 629 (1990) 

39. D.P. Barber, G. Ripken, DESY HERA Report 87-15 (1987) 

40. J. Buon, Proc. 6th Int. Symp. High Energy Spin Physics, Journal de Physique, 
Colloque C2, Tome 46 (1985) 

41. D.P. Barber, Private notes (1995) 

42. K. Yokoya, Particle Accelerators 13 (1983) 85 

43. S.R. Mane, Nucl.Inst.Meth. A292 (1990) 52 



22 



44. S.R. Mane, Nucl.Inst.Meth. A321 (1992) 21 

45. Handbook of Accelerator Physics and Engineering, Eds. A.W. Chao and M. 
Tigner, World Scientific (1999) 



23 



